Doped carrier formulation and mean-field theory of the tt't" J model 
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In the generalized-t J model the effect of the large local Coulomb repulsion is accounted for by 
restricting the Hilbert space to states with at most one electron per site. In this case the electronic 
system can be viewed in terms of holes hopping in a lattice of correlated spins, where holes are the 
carriers doped into the half-filled Mott insulator. To explicitly capture the interplay between the 
hole dynamics and local spin correlations we derive a new formulation of the generalized-f J model 
where doped carrier operators are used instead of the original electron operators. This "doped 
carrier" formulation provides a new starting point to address doped spin systems and we use it to 
develop a new, fully fermionic, mean-field description of doped Mott insulators. This mean-field 
approach reveals a new mechanism for superconductivity, namely spinon-dopon mixing, and we 
apply it to the tt't" J model as of interest to high-temperature superconductors. In particular, we 
use model parameters borrowed from band calculations and from fitting ARPES data to obtain 
a mean-field phase diagram that reproduces semi-quantitatively that of hole and electron doped 
cuprates. The mean-field approach hereby presented accounts for the local antiferromagnetic and 
d-wave superconducting correlations which, we show, provide a rational for the role of t' and t" 
in strengthening superconductivity as expected by experiments and other theoretical approaches. 
As we discuss how t, t' and t" affect the phase diagram, we also comment on possible scenarios to 
understand the differences between as-grown and oxygen reduced electron doped samples. 
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I. INTRODUCTION 

High-temperature superconducting (SC) cuprates are 
layered materials where the copper-oxide planes are sep- 
arated by several elements whose chemistry controls the 
density of carriers in the Cu02 layers. This density deter- 
mines the electronic nature of copper-oxide planes, where 
the renowned unconventional cuprate phenomenology is 
believed to take place. Due to the localized charac- 
ter of 3d orbitals, copper valence electrons feel a large 
Coulomb interaction which drives strong correlations .A 
The largest in-plane energy scale is the local Coulomb 
repulsion U, measured to be approximately 2eVjfi thus 
motivating the use of the generalized-^ J model in the 
cuprate context Interestingly, upon inclusion of pa- 
rameters t' and t" consistent with electronic structure 
calculations^ the generalized-t J model reproduces many 
spectral features observed in both hole and electron 
doped cupratesAi 

Following its aforementioned relevance, in this paper 
we consider the two-dimensional tt't" J Hamiltonian 



H 



tJ 



j £ (s, Sj -\ 
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where £y = t, t' , t" for first, second and third nearest 
neighbor (NN) sites respectively. c\ 



- [ct |ct , ] is the elec- 
: " and Sj = cltrcj 



tron creation spinor operator, rij = c\Ci 
are the electron number and spin operators and er are 
the Pauli matrices. The operator V projects out states 



with on-site double electron occupancy and, therefore, 
the tt't" J model Hilbert space consists of states where 
every site has either a spin-1/2 or a vacancy. Conse- 
quently, when the number of electrons equals the num- 
ber of lattice sites, i.e. at half-filling, there exists ex- 
actly one electron on each site and the system is a Mott 
insulator which can be described in terms of spin vari- 
ables alone. In the doped Mott insulator case, how- 
ever, the model Hamiltonian also describes electron hop- 
ping, which is highly constrained since electrons can only 
move onto a vacant site. This fact is captured by the 
use of projected electron operators Vc? irj V that do not 
obey the canonical fermionic anti-commutation relations 
and, thus, constitute a major hurdle to handle the tt't" J 
model analytically. For that reason, the physics of the 
generalized-t J model has been addressed by a variety of 
numerical techniques, which include the exact diagonal- 
ization of small systemsjiSuLifiiiiii^iiii the self-consistent 
Born approximationji^iSiiSiiLiSiiS the cellular dynami- 
cal mean-field (MF) approximation^ quantum Monte 
Carlo, 9 i 21 i 22 i 23 i 24 Green function Monte Carlo 2 ^ and 
variational Monte Carlo methodsi 27 i 28 i 29 i 30 These stud- 
ies provide an overall consistent picture that guides the 
effort to develop analytical approaches to the generalized- 
t J model, which is the main focus of the current work, 
and we refer to them throughout the paper. 

Different analytical techniques and approximations 
have been developed to address doped Mott insulators in 
terms of bosonic and/or fermionic operators amenable to 
a large repertoire of many-body physics tools. The sim- 
plest one, known as the Gutzwiller approximationj 3 ^^ 
trades the projection operators in Vc\cjV by a nu- 
merical renormalization factor g t — that vanishes 
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in the half-filling limit. The extensively used slave- 
particle techniqueo 14 ! 33 ! 34 . 35 ! 36 ! 37 ! 38 ! 39 ! 40 ! 41 rather decou- 
ple the projected electron operator Vc\ a V into a fermion 
and a boson which describe chargeless spin- 1/2 excita- 
tions of the spin background and spinless charge e exci- 
tations that keep track of the vacant sites. These various 
approaches provide distinct ways to handle the projected 
electron operator Vc\ a V which, we remark, significantly 

differs from the bare electron operator c\ close to half- 
fillingi 2 . 

In this paper, we present and explicitly derive a differ- 
ent approach to doped Mott insulators that circumvents 
the use of projected electron operators. This approach 
was first introduced in Ref. [43| and recasts the tt't" J 
model Hamiltonian Q in terms of spin variables and pro- 
jected doped hole operators which carry the charge and 
spin degrees of freedom introduced in the system upon 
doping. To make the difference between both approaches 
more concrete, note that while projected electron opera- 
tors are used to fill the sea of interacting electrons start- 
ing from the empty vacuum, projected doped hole oper- 
ators are used to describe the sea of interacting carriers 
doped into the half-filled spin system. 

One should bear in mind that the physical properties of 
holes in doped Mott insulators differ from those of holes 
in conventional band insulators. The distinction is clear 
in the infinite on-site Coulomb repulsion limit, i.e. in the 
generalized-^ J model regime, in which case the insertion 
of a hole on a certain site i amounts to the removal of 
a lattice spin. This leaves a vacant site, that carries a 
unit charge e when compared to the remaining sites, and 
changes the total spin component along the z-direction 
by ±^4^ Since the vacancy is a spin singlet entity, the ex- 
tra spin ± i introduced upon doping is not on the vacant 
site and is carried by the surrounding spin background 
instead. Hence, in doped Mott insulators, the object 
that carries the same quantum numbers as doped holes, 
namely charge e and spin ±i, is not an on-site entity^ 
We remark that in the doped Mott insulator literature 
the term "hole" is often used with different meanings. 
Sometimes it refers to the charge-e spin-0 object on the 
vacant site from which a spin was removed. We reserve 
the term "vacancy" for such an on-site object which dif- 
fers from the above described entity which carries both 
non-zero charge and spin quantum numbers. We use the 
term "doped hole" to allude to this latter non-local en- 
tity. In fact, and since the tt't" J model can be used to 
describe both the lower Hubbard band or the upper Hub- 
bard band in the large Coulomb repulsion limit, in this 
paper we generally use the term "doped carrier" to mean 
"doped hole" or "doped electron" depending on whether 
we refer to the hole or electron doped case. 

It follows from the above argument that a doped car- 
rier in doped Mott insulators is a composite object that 
involves both the vacancy and its surrounding spins 
and, therefore, its properties are related to those of the 
background spin correlations. Since at half-filling the 



generalized-i J model reduces to the Heisenberg Hamil- 
tonian, which displays antifcrromagnetic (AF) order, 46 
it is natural to think of doped carriers in terms of a va- 
cancy encircled by a staggered spin configuration. As 
emphasized in Ref. '261, this picture is valid even in the 
absence of long-range AF order and, in fact, it only re- 
quires that most of the doped carrier spin-1/2 is concen- 
trated in a region whose linear size is smaller than the AF 
correlation length. Interestingly, quantum Monte Carlo 
calculations for the U/t = 8 Hubbard model^i find signa- 
tures of short-range AF correlations around the vacancy 
up to the hole density x — 0.25. Angle-resolved pho- 
toemission spectroscopy (ARPES) experiments also show 
that the high energy hump present in undoped samples, 
which disperses in accordance with the tt't" J model sin- 
gle hole dispersion;£ii2iiiiiiii is present all the way into 
the overdoped regime4£ Hence, both theory and exper- 
iments support that short-range AF correlations persist 
throughout a vast range of the high-T c phase diagram 
and, thus, suggest that the above local picture of a doped 
carrier holds as we move away from half-filling. 

It is important to understand how superconductivity 
arises out of the above doped carrier objects. One pos- 
sibility is that doped carriers form a gas of fermionic 
quasiparticlesi 3 ^ which form local bound states in the 
d-wave channel due to AF correlations. 26 Alternatively, it 
has been proposed- that doped carriers induce spin liquid 
correlations which do not frustrate the hopping of charge 
carriers (unlike AF correlations) and which ultimately 
lead to SC order. ARPES data on the cuprates identi- 
fies two different dispersions - at low energy the char- 
acteristic d-wave SC linear dispersion crosses the Fermi 
level close to (-§,§) while at higher energy the disper- 
sion inherited from the undoped limit persists 
This experimental evidence supports the coexistence of 
both local AF and d-wave spin singlet correlations and, 
thus, favors the second scenario above. The coexistence 
of two distinct types of spin correlations also receives sup- 
port from exact diagonalization calculations of the tt't" J 
models 

Since the formulation of the tt't" J model hereby pre- 
sented introduces doped carrier operators to address 
doped Mott insulators close to half-filling we dub it the 
"doped carrier" formulation of the tt't" J model. A num- 
ber of advantages stem from using doped carrier rather 
than the original particle operators. For instance, close 
to half-filling the doped carrier density is small and so 
we are free of the no-double-occupancy constraint prob- 
lem for doped carriers. In addition, as we show in this 
paper, in the "doped carrier" framework the hopping sec- 
tor of the tt't" J model explicitly describes the interplay 
between the doped carrier dynamics and different local 
spin correlations, specifically the aforementioned coexist- 
ing local AF and d-wave SC correlations. Consequently, 
this approach provides a powerful framework to under- 
stand the single-particle dyna mics in doped Mott insu- 
lators, as attested by Refs. I43ll52l53l which show that it 
can be used to reproduce a large spectrum of ARPES 
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and tunneling spectroscopy data on the cuprates. 



In Sec. |n] we present the detailed derivation of the 
"doped carrier" formulation of the tt't" J model. In par- 
ticular, we define the aforementioned projected doped 
carrier operators, which we call "dopon" operators, as 
on-site fermionic operators with spin-1/2 and unit elec- 
tric charge. Dopons act on an Hilbert space larger than 
the tt't" J model Hilbert space and thus they are distinct 
from electron operators. Only after the contribution from 
all unphysical states is left out do dopons describe the 
physical doped carriers, which are the non-local entities 
consisting of the vacancy and the extra spin ± ^ carried 
by the encircling spins. 



The "doped carrier" framework constitutes a new 
starting point to address doped spin systems and in Sec. 
IIIII we construct a new, fully fermionic, MF theory of 
doped Mott insulators, which we call the "doped car- 
rier" MF theory. In this MF approach we represent 
the lattice spin variables, which provide the half-filled 
background on top of which doped carriers are added, in 
terms of fermionic spin-1/2 chargeless excitations known 
as spinons. Hence, we describe the low energy dynamics 
of the tt't" J model in terms of spinons and dopons. Since 
the latter carry the same quantum numbers as electrons, 
the "doped carrier" formulation can account for low ly- 
ing electron-like quasiparticles which form Fermi arcs 
as observed in high temperature superconductors. 43 - 52 - 54 
This state of affairs is to be contrasted with that stem- 
ming from the slave-boson approach to the generalized-t J 
modeli 33 i 34 i 37 i 38 i 39 i 40 i 41 In the latter approach the low en- 
ergy dynamics is described by spinons and holons (spin- 
less charge-e bosons) and, thus, it stresses the physics of 
spin-charge separation. Despite the differences, it is pos- 
sible to relate the "doped carrier" and the slave-boson 
formulations of the generalized-f J model and, in Ap- 
pendix we show that holons are the singlet bound 
states of spinons and dopons. 



In Sec. II VI we use model parameters extracted both 
from ARPES data and from electronic structure cal- 
culations concerning the copper-oxide layers of high- 
temperature superconductors to discuss the MF phase 
diagrams of interest to hole and electron doped cuprates. 
The relevance of the hereby proposed approach is sup- 
ported by the semi-quantitative agreement between our 
results and the experimental phase diagram for this fam- 
ily of materials. The "doped carrier" MF theory explic- 
itly accounts for how local AF and <i-wave SC correlations 
correlations affect and are affected by the local doped car- 
rier dynamics and, thus, we are able to discuss the role 
played by the hopping parameters t, t' and t" in deter- 
mining the MF phase diagram. These and other results 
underlie our conclusions (Sec. Ivjl . 



II. DOPED CARRIER FORMULATION OF 
THE tt't" J MODEL 

A. Enlarged Hilbert space 

The tt't" J model Hamiltonian is written in terms 
of projected electron operators Pc] a V and Vci^V which 
rule out doubly occupied sites and, thus, the tt't" J model 
on-site Hilbert space for any site i is 

Wi = {|T) i ,|i) i ,|0) i } (2) 

The states in Q include the spin-up, spin-down and va- 
cancy states respectively. 

In this section, we introduce a different, though equiva- 
lent, formulation of the tt't" J model. In the usual formu- 
lation any wave-function can be written by acting with 
the projected electron creation operators Vc\ a V on top 
of an empty background. Close to half-filling these oper- 
ators substantially differ from the bare electron creation 
operator and we propose an alternative framework where 
the background on top of which carriers are created is 
a lattice of spin-1/2 objects, which we call the "lattice 
spins". We then consider fermionic spin-1/2 objects with 
unit charge, which we call "dopons", that move on top 
of the lattice spin background. Dopons have the same 
spin and electric charge quantum numbers as the carri- 
ers doped in the system and are introduced to describe 
these doped carriers. In such a description there is one 
lattice spin in every site whether or not this site corre- 
sponds to a physical vacancy. In addition, there exists 
one, and only one, dopon in every physically vacant site. 
However, a vacant site is a spinless object while dopons 
carry spin-1/2. Therefore, to accommodate the presence 
of both lattice spins and dopons we must consider an 
enlarged on-site Hilbert space which for any site i is 

w? B, = {|T0) i ,U0) i ,|U) iJ UT) iJ |Tt) i ,|U) i } (3) 

States in © are denoted by \aa), where a =|,J. labels 
the up and down states of lattice spins and a = 0, T, J. 
labels the three dopon on-site states, namely the no do- 
pon, spin-up dopon and spin-down dopon states. To 
act on these states we introduce the lattice spin opera- 
tor Si and the fermionic dopon creation spinor operator 
d\ — [dJi-dJjJ, which are such that Sf \aa) i = s a \ \aa) i 
[where s a = (+) for a =| and s a = (— ) for a =]\ 
and d\ a \aO) i = \aa). Since in the enlarged Hilbert 
space (0 there exist no states with two dopons on 
the same site we also introduce the projection operator 
V = rii(l — d\ i-di,f<4 j^i.i) which enforces the no-double- 
occupancy constraint for dopons. 

In order to write physical operators, like the projected 
electron operators or the tt"t" J model Hamiltonian 0J, 
in terms of lattice spin operators Si and projected do- 
pon operators Vd\ a V and Vdi^ a V we define the follow- 
ing mapping from states in the enlarged on-site Hilbert 
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space (0 onto the physical tt't" J model on-site Hilbert 
space © 



It); - ITO), ; \ih - 140), ; 10), 



IU)i-|lt)i 

IU>»+I4T>< 



(4) 



The on-site triplet states in Hf n , namely nih^dli^ |||) . 

and III), are unphysical as they do not map onto any 
state pertaining to the tt't" J model on-site Hilbert space. 
Therefore, these latter states must be left out when writ- 
ing down the physical wave-function or when defining 
how physical operators act on the enlarged Hilbert space. 



B. Electron operator in enlarged Hilbert space 



The mapping rules (@J) can be used to define the spin- 
up electron creation operator c] ^ that acts on the on- 
site enlarged Hilbert space and whose matrix elements 
on the physical sector of J3J) match those of Vc\ on 
(J5J), namely 



v<{<(P{W)i\\l)i\m = {0;0;|T)J 



(5) 



If we further require cj ^ to vanish when it acts on un- 
physical states we obtain the relations 



-t [Ti>« - UT>< 



V2 



= Ito), 



,t ITO), = < T IIO) i; - 
4 T (Itl)i + lit);) = < T ltt)i = e], T IU>< = o (6) 

It is a trivial matter to recognize that the defining rela- 
tions lO are satisfied by the operator 



C i,T 



V2 



S+d 



(7) 



where fif = Sf±iSf. 

The spin-down electron creation operator can be de- 
termined along the same lines or, alternatively, by con- 
sidering an angle 7r rotation around the y-axis, which 
transforms 



5 f 

Sf 



-5 f 



HA 



-S? ; S; 



~ d i,l 

-sr 



sr - 



di,t 

-s+ 



(8) 



It then follows that the electron operatorsSS that act 
on the on-site enlarged Hilbert space are 



5 - = ^ 



1 



Sf'ek. 



(9) 



C. tt't" J model Hamiltonian in enlarged Hilbert 
space 

We now recast the tt't" J model Hamiltonian QJ, which 
is a function of the projected electron operators "Pcj a V 



and Vci !<T V, in terms of lattice spin and projected dopon 
operators as dictated by Expression ©. The resulting 
tt't" J model Hamiltonian in the enlarged Hilbert space 



Henl ttJ . 77 
tj — -rc„»,7 + tl. 



J , TTt 

enl ' enl 



(10) 



is the sum of the Heisenberg {H^ nl ) and the hopping 
{Herd) terms. First we use the equalities 

Vd^Vd\_J> = V{d iA Vd\ A - di, L Vdl d )P = 
V{d lA Vd\ A + d lA Vd\ A )V = 2V(l - d\di)V (11) 

to replace the operators Si = {Vc i i 'P)(T{'Pc{P) and 
VriiV = (T4v)(VciT) in the first term of Q by 

S{P{1 - d\di)V and V(l - d\<k)V respectively. As a 
result, the Heisenberg interaction in terms of lattice spin 
and projected dopon operators is 



J £ U.S., \ | 



(l - d\di) (l 



I - d\dj ) T 



(12) 

We obtain the hopping term in the enlarged Hilbert space 
upon directly replacing Vc\ a V and Vc^ a V in the second 
term of QJ by c\ and Cj )(T , which leads to 



(ij) 




1 S i X S 7 _ 



1 



djdj + dld.jSi.Sj + h.c. 



P (13) 



The hopping term in the tt't" J model Hamiltonian 
connects a state with a vacancy on site j and a spin 
on site i to that with an equal spin state on site j and a 
vacancy on site i. This is a two site process that leaves the 
remaining sites unaltered and, schematically, it amounts 
to 



|0;T) 
|0;l) 



IT; o> 
II; o) 



(14) 



where the notation \j; i) is used to represent the states on 
sites j and i. Using the corresponding two site notation 
for the enlarged Hilbert space, namely \ajCij] o^a;), where 
cti,ctj =T,I and a.i,aj = 0, |,|, and making use of the 
mapping rules |4f. the hopping processes in the "doped 
carrier" framework that correspond to (|14|l are 



Tl - IT 



V2 

TI-IT 



V2 



;|0 

;|o 



|0: 
10: 



TI-IT 
TI-IT 

V2 



(15) 



It can be explicitly shown that the hopping term 113|) 
is such that its only non-vanishing matrix elements are 
those that describe the processes in (|15H . Hence, only 
local singlet states hop between different lattice sites 
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whereas the unphysical local triplet states are localized 
and have no kinetic energy. Therefore, the dynamics de- 
scribed by H l enl effectively implements the local singlet 
constraint, which leaves out the unphysical states in the 
enlarged Hilbert space ©. 

We emphasize that the tt't" J Hamiltonian in the en- 
larged Hilbert space l|l(J equals H t j in the physical 
Hilbert space. In addition, it does not connect the phys- 
ical and the unphysical sectors of the enlarged Hilbert 
space. Therefore, the "doped carrier" formulation of the 
tt't" J model, as defined by (fTTTf) - fT2*)l and fT3|) . is equiva- 
lent to the original "particle" formulation encoded in Q . 
Interestingly, it provides a different starting point to deal 
with doped spin models. 

In the low doping regime the dopon density x — 
jf^iVdldiV, where N is the number of sites, is small 
and the no-double-occupancy constraint for dopons is 
safely relaxed. Hence, below we drop all the projection 
operators V. We thus propose that the dramatic effect of 
the projection operators V in the "particle" formulation 
of the tt't" J model Q is captured by the dopon-spin in- 
teraction in the hopping hamiltonian (|13|l . which explic- 
itly accounts for the role of local spin correlations on the 
hole dynamics. In the remainder of the paper we derive 
and discuss a MF theory that describes this interaction. 



III. DOPED CARRIER MEAN-FIELD THEORY 
OF THE tt't" J MODEL 

To derive a MF theory of the tt't" J model "doped car- 
rier" formulation we recur to the fermionic representa- 

is 



tion of lattice spins Si = |//tr/i, where fj = [/^/^j 
the spinon creation spinor operator^ The Hamiltonian 
Hlf is then the sum of terms with multiple fermionic 
operators which can be decoupled upon the introduc- 
tion of appropriate fermionic averages, as presented in 
what follows. The resulting MF Hamiltonian H^j F is 
quadratic in the operators f<, /, d^ and d, and de- 
scribes the hopping, pairing and mixing of spinons and 
dopons. We remark that, in contrast to slave-particle ap- 
proaches which split the electron operator into a bosonic 
and a fermionic excitationsji^i&iil this MF Hamiltonian 
is purely fermionic. 



A. Heisenberg term 

We first consider the spin exchange interaction i|12|) 
which upon replacing the operator V(l — d\di)V by its 
average value (1 — x) reduces to the Heisenberg Hamil- 
tonian 



(16) 

(ij)eNN 

with the renormalized exchange constant J — (1 — x) 2 J. 



In the enlarged Hilbert space there is always one 
lattice spin Si per site, even in the presence of finite 
doping. Therefore, in the fermionic representation of Si 
the SU(2) projection constraint enforcing f\fi = 1 must 
be implemented by adding the term eio.i- (Vl <T V'i) to 

(|10|l . where a^ : i are local Lagrangian multipliersi 3 ^ 4 ^ 
Here we use the Nambu notation for spinon operators, 
namely ip\ = [V'liV^] = L/i j/«,lL m terms of which 
lattice spin operators can be recast as Si.cr = K^/tv]^— I) 
where 3 ^ 



V, = 



4,2 

-4a 



(17) 



and I is the identity matrix. It then follows that, in 
the fermionic representation, the Hamiltonian term 
reduces to 



J 



(ij)(=NN 



cto,, 



(18) 



with U id = -*j*t. 

The quartic fermionic terms in l|18|l can be decoupled 
by means of the Hartree-Fock-Bogoliubov approximation 
thus leading to the MF Heisenberg terr 1 33 ' 34 



H 



MF 



3J 
16 



3 J 



Tr [Ui tj U jti ] +a .lj2 ^Wi ) 

)GNN V i J 



(19) 



(ij)GNN 



where Ui.j — (Ujj) are the singlet bond MF order param- 
eters (in Sec. IIII Fl the above MF decoupling is extended 
to include the formation of local staggered moments as 
well). Also note that at the MF level the local con- 
straint is relaxed and ao is taken to be site independent. 
The best non-symmetry breaking MF parameters to de- 
scribe the paramagnetic spin liquid state of the Heisen- 
berg model correspond to the ansatz 

Ui.i+x = XOz + Acr x ; Ui^i+y = \a z - Aa x ; a = a Q a z 

(20) 

which describes spinons paired in the d-wave channel and 
whose properties have been studied in the context of the 
slave-boson approachi 3 ^ 3 ^ 4 * 41 . 



B. Hopping term 

We now consider the hopping Hamiltonian l|13l) , which 
describes the hopping of holes on the top of a spin back- 
ground with strong local AF correlations. It is well un- 
derstood that, under such circumstances, the hole disper- 
sion is renormalized by spin fluctuations Ai^S^S As it is 
further elaborated in Sec. IIII El to capture the effect of 
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this renormalization at the MF level we replace the bare 
t, t' and t" by effective hopping parameters t\, t 2 and i 3 . 

The hopping term (fT3*|l involves both spinon operators 
(1 1 71) and dopon operators 



Di 



a i,l 

-4t 



(21) 



and after dropping the projection operators V it can be 
recast as 

*f [tt [(*}*.,- _ /) Dj<j z Di - I 



-Tr 
-Tr 



D]a z Di 



Tr 



(22) 



The contribution from l|22|) to the MF Hamiltonian, 
namely H\^ F , can be decomposed into the MF terms 
that belong to the spinon sector (H t M s F nnon ), the ones 
that contribute to the dopon sector (H t £j F \ pon ) and the 
ones that mix spinons and dopons (H^p*), so that 

Trt,dopon rrt,mix „ , r 

a if p i -tJ a i j? • rui liii^ oai\.c ui 



n MF 



jjt. spinon 
n MF 



MF t ±± MF 

clarity, in what follows, we discuss each of the above three 
contributions separately: 

(i) H M S p mon arises from decoupling the first term in 
(J22J- The decoupling is done by taking the average 

(^Djtr.AtfJ) - (B)^ z B l;l ) « xa z 6 meNN (23) 



where we introduce the spinon-dopon singlet pair opera- 



tor /;,, = Dj^\. Then 60 



rj-t, spinon 
MF 



E f(4^+h. 



(24) 



(ij)ENN 



This term determines the effect of spinon-dopon pairs 
Bj^i in the magnitude of the spinon d-wave pairing gap. 
As we discuss in Appendix ^ the spinon-dopon sin- 
glet pair operator Bji corresponds to the holon oper- 
ator in the slave-boson formalism and, in fact, a term 
similar to l|24[) appears in the SU(2) slave-boson MF 
Hamiltonian^i 

(ii) The MF dopon hopping term i?^ po ™ comes from 
the fourth term in (|22|l and from taking the average 



*J*,-Z ))=A(S i .S j +icr. (SixSj ) 



(25) 



in the first term of l|22|) . As we mention in Sec. [I] both 
numerical 9 i 13 i 21 i 24 i 25 i 26 and experimental^ 48 ! 49 ! 50 ! 51 evi- 
dence point to the persistence of signatures due to local 
AF correlations around the vacancies all the way into the 
doping regime where cuprates superconduct. To account 
for this effect we consider that the spins encircling the 
vacancy in the one-dopon state are in a local Neel con- 
figuration and, therefore, in l|25|) we use (Si x Sj) =0 



and A(Si.Sj) = (— As a result, the contri- 
bution from the first and fourth terms in (1221) to the MF 
dopon sector is 



H 



t : dopon 

MF 



= E t (4°*^ + h 



(ij)£2 nd NN 



E jfa^ + h.c) (26) 



{ij)£3 rd NN 



where the dopon Nambu operators m = TO 

are introduced. Ideally, the average in f2*5|) 
should be calculated self-consistently to reproduce the 
doping induced change in the local spin correlations. In 
the present MF scheme this is not performed and, in- 
stead, we introduce doping dependent effective hopping 
parameters t 2 and t 3 (see Sec. MI E|) . 

(in) The MF term that mixes dopons and spinons, 
^mf %X i captures the interaction between the spin degrees 
of freedom and the doped carriers enclosed in the second 
and third terms of l|22|) . which can be recast as 



t l = 



E 



-Tr 



(27) 



In the Hartree-Fock-Bogoliubov approximation the above 
expression leads to the spinon-dopon mixing term 

H'mT = - E Tr - E [4Bii^i + h.c 



16 



E*" E U+uBi^i+h.c) (28) 



uBv NN 



where B;, 



B iti and B lX 



16 Tlv Efiew NN Bi+u^i- 

Here, the mean-fields Bj A = (BaJ) are introduced. We 
also use u = ±x, ±y, u — ±x ± y and u = ±2x, ±2y for 
v = 1,2,3 respectively. Different choices for the mean- 
fields B i0 and Bn may describe distinct MF phases. In 
what follows we take 



Bm = —boo- z 



Bi 



-bxa z 



(29) 



which, as we show in Appendix [S] describes the d-wave 
SC phase when taken together with the spinon <i-wave 
ansatz (j20J). 

In order to clarify the physical picture enclosed in 
the above MF scheme, note that dopons correspond 
to vacancies surrounded by a staggered spin configu- 
ration and, therefore, they describe quasiparticles in 
the half-filling limit. Such a locally AF spin back- 
ground strongly suppresses coherent doped carrier inter- 
sublattice hopping 4 *^ 5 ^ 9 . as captured by H M d p P ° n which 
includes dopon hopping processes between 2 nd and 3 rd 
NN sites but not between 1 st NN sites. At MF level, 
the only mechanism for dopons to hop between differ- 
ent sublattices is provided by the spinon-dopon mixing 
term l|28|) which represents the interaction between do- 
pons and the lattice spins. This interaction leads to the 
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formation of spinon-dopon pairs -Bj.i, which are spin sin- 
glet electrically charged objects and, thus, describe va- 
cancies surrounded by spin singlet correlations that en- 
hance the hopping of charge carriers. From l|29[) we have 
that b = (fjdi) and h = £„ £ fiG „ NN /M+fi) 
are the local and non-local MF parameters that emerge 
from such spin assisted doped carrier hopping events and 
which describe the hybridization of spinons and dopons. 
In the rest of the paper we interchangeably refer to the 
condensation of the bosonic mean-fields bo and b\ as the 
coherent spinon-dopon hybridization, mixing or pairing. 
As a final remark, note that for bij^O the term ~ bid] ft in 
(|28J) drives local mixing of spinons and dopons and leads 
to a non-zero bo- Similarly, if bo^O the term ~ bod\,^fi 
leads to non-local spinon-dopon mixing and, thus, to non- 
zero b\. Hence, either bo and b\ are both zero or both 
non-zero. 



Mean-field Hamiltonian 



Putting the terms (JUJ), JUJ), (J2SJ and together 
leads to the full "doped carrier" MF Hamiltonian H^ F 



H J 

MF 



HJ 



which in momentum space becomes 



a%a z + a%a x (3 k a z 

Pk&z lk&z 





4>k 




Vk 



3JN 



(X 2 + A 2 ) - 2Nb h - Nn d 



(30) 



where 



3J 



X — t\x (cos k x + cos k y ) + ao 



3J A / , 

A (cos k x 

4 ^ x 



3bo 



/3fe = — [h (cosk x + cos fc y ) 
o 



COS fe, 



2to COS fer COS fe, 



7k = t2 cos k x cos ky + — (cos 2k x + cos 2k y ) — fid (31) 



+t 3 (cos 2k x + cos 2k y )] 

h 

2 

In (|30|l we introduce the dopon chemical potential fid 
that sets the doping level (d\di) = x. The explicit form 
of the above MF Hamiltonian depends on the values of 
ti, £2 and £3, which are determined phenomenologically 
in Section fill El by fitting to both numerical results and 
cuprate ARPES data. The mean-field parameters x, A, 
bo, and b\ are determined by minimizing the mean-field 
free-energy and in Section IIVI we show they reproduce 
the cuprate phase diagram. 

D. Two-band description of one-band tt't" J model 

Even though the tt't" J model is intrinsically a one- 
band model, the above MF approach contains two dif- 
ferent families of spin- 1/2 fermions, namely spinons and 



dopons, and thus presents a two-band description of the 
same model. As a result, H^j has a total of four 
fermionic bands described by the eigenenergies 



c i,fc 



-2,fc 



±ypk 



±ypk 



(32) 



where 



1 

2 



it K) - Kr 



Pk 



it + K) + Kr 



(33) 



In the absence of spinon-dopon mixing, i.e. when 
60,61 = 0, the bands e^ fc and e 2 k describe, on the one 
hand, the spinon <i-wave dispersion that underlies the 
same spin dynamics as obtained by slave-boson theorypS 
In addition, these bands also capture the dispersion of 
a hole surrounded by staggered local moments which in- 
cludes only intra-sublattice hopping processes [Expres- 
sion (|2*oTl ]. as appropriate in the one- hole limit of the 
tt't" J modelA 2 ^^ Upon the hybridization of spinons 
and dopons the eigenbands ef k and e^ k differ from 
the bare spinon and dopon bands by a term of or- 
der 6q, 61,60^1 ~ x. In particular, the lowest energy 
bands t x k are d-wave-like with nodal points along the 
(0, 0) — (±7r, ±7r) directions and describe electronic exci- 
tations that coherently hop between NN sites. The high- 
est energy bands e 2 k are mostly derived from the bare 
dopon bands and, therefore, describe excitations with re- 
duced NN hoppingi 43 ! 52 

The reason underlying the above multi-band descrip- 
tion of the interplay between spin and local charge dy- 
namics stems from the strongly correlated nature of the 
problem. Quasiparticles in conventional uncorrelated 
materials correspond to dressed electrons whose disper- 
sion depends to a small extent on the remaining excita- 
tions and is largely determined by an effective external 
potential. In the presence of strong electron interactions, 
though, the dynamics of electronic excitations is inti- 
mately connected to the surrounding environment and 
depends on the various local spin correlations. Physically, 
the two-band description provided by the "doped carrier" 
MF theory captures the role played by two such differ- 
ent local spin correlations on the hole dynamics. These 
are the local staggered moment correlations, which are 
driven by the exchange interaction but frustrate NN hole 
hopping, and the d-w&ve spin liquid correlations, which 
enhance NN hole hopping at the cost of spin exchange 
energy. 

We remark that such a multi-band structure agrees 
with quantum Monte Carlo and cellular dynamical MF 
theory calculations on the two-dimensional t J and Hub- 
bard models. 9 i 20 i 21 i 22 i 23 In Refs. I22l23l the two bands 
below the Fermi level were interpreted in terms of two 
different states, namely: (i) holes on the top of an other- 
wise unperturbed spin background and (ii) holes dressed 
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by spin excitations. This interpretation is consistent with 
additional numerical work indicating the existence of two 
relevant spin configurations around the vacancj^ and of- 
fers support to the above MF formulation. 

Finally, we point out that the two-band description 
of the generalized-i J model "doped carrier" formulation 
resembles that of heavy-fermion models: dopons and lat- 
tice spins in the "doped carrier" framework correspond 
to conduction electrons and to the spins of /-electrons, 
respectively, in heavy-fermion systems. The main dif- 
ference is that, at low dopings, the tt't" J model spin- 
spin interaction is larger than the dopon Fermi energy, 
while in heavy-fermion models the spin-spin interaction 
between /-electrons is much smaller than the Fermi en- 
ergy of conduction electrons. As the doping concentra- 
tion increases the dopon Fermi energy approaches, and 
can even overcome, the interaction energy between lat- 
tice spins. In that case, our approach to the tt't" J model 
becomes qualitatively similar to heavy-fermion models. 
Interestingly, overdoped cuprate samples do behave like 
heavy-fermion systems except for the mass enhancement, 
which is not as large as for typical heavy-fermion com- 
pounds. 



E. Renormalized hopping parameters 

In Sec. IIII Bl we mention that at the MF level we resort 
to effective hopping parameters t\, t 2 and £3 to account 
for the renormalization due to spin fluctuations. In this 
paper we take the NN hopping parameter to equal its 
bare value, i.e. t\ = t. However, as we discuss in what 
follows, the role of local spin correlations on the intra- 
sublattice doped carrier dynamics is quite non-trivial and 
t 2 and £3 differ from the corresponding bare parameters. 



1. Single hole limit 

The "doped carrier" MF theory of the tt't" J model 
considers the dilute vacancy limit where we can focus 
on the local problem of a single vacancy surrounded by 
spins. In particular, it captures the effect of strong 
AF correlations around the vacancy through the av- 
erage H25J) which determines the dopon dispersion 7^. 
This dispersion is controlled by the hopping parame- 
ters t 2 and t 3 and does not involve NN hopping pro- 
cesses, in agreement with tt't" J model single hole prob- 
lem results l ^ 1 9i l; ?' 1 ^ 1 Ti 1 ?i 2; Pi 2 ?i^ However, the values of 
t 2 and £3 that fit the above single hole dispersion dif- 
fer from the bare t' and t" . This is clearly so when 
t' = t" = and, still, charge carriers move coherently 
within the same sublatticeiiiiS* 2 ^ These intra-sublattice 
hopping processes, which induce non-zero effective hop- 
ping parameters t 2 and t%, result from a spin fluctua- 
tion induced contribution to charge hopping^ We re- 
mark that this is inherently a quantum contribution that 
escapes the realm of MF theory, which is a semiclassical 



saddle-point approach, and therefore we recur to numer- 
ical and experimental evidence to set the values of t 2 and 

It is well established that for t! — t" = the single hole 
dispersion has its minimum at {tt/2, tt/2) and is quite flat 
along (0, tt) — (tt, 0) pfiiifliiLiSiS 5 . For our purposes, we can 
simply consider that the hole dispersion is completely 
flat along (0, 7r) — (tt, 0) and, thus, we set t 2 = 2t 3 = tj 
when t' and t" vanish. In this case, numerical calcula- 
tions show that the dispersion width along (0, 0) — (tt, tt) 
is w 2J 4 -' 6 and we take tj = J. The resulting MF effective 
hopping parameters t 2 — J and £3 = 0.5 J compare well 
with those found by the self-consistent Born approxima- 
tion for t = 0.3J, namely t 2 = 0.87J and t 3 = 0.62J^£ 
and with those determined by the Green function Monte 
Carlo technique for t = 0.4J, specifically t 2 = 0.85J and 
t 3 = 0.65J-2 5 . 

Since non-zero values of t' and t" do not frustrate nor 
are frustrated by local AF correlations, in this case we 
naively take the effective intra-sublattice hopping param- 
eters introduced in Sec. IIII Bl to be 

t 2 = tj + t' 

^ = Y+t" (34) 

Based on experimentsi2*§£ and band theory calculations 5 
relevant to the cuprates which show that bare intra- 
sublattice hopping parameters are non-zero together 
with numerical calculation results that support t' ~ 
_ 2^6,10,15,16 we set fj = 2t". The effective hopping pa- 
rameter choice in (|34H then leads to a dopon dispersion 
width along (0, 7r) — (71", 0) equal to 2t' , which given the 
simple approximation scheme involved is reasonably close 
to numerical results findings, namely that the dispersion 
width along (0, tt) — (tt, 0) is At' where the coefficient A 
is somewhere in the range | — 4^ 

We emphasize that the above comments to the sin- 
gle hole dispersion in the tt't" J model appear to be rel- 
evant to the cuprate materials. Indeed, a large body 
of experimental evidence shows that the nodal disper- 
sion width as measured by ARPES in undoped samples 
is ft 2J) 1 9i 4 7i 4 9i 4 ?iWi 6i i l where J is independently deter- 
mined by band calculations^ and from fitting Raman 
scattering^ and neutron scattering^ experiments. In 
addition, the experimental dispersion along (0, tt) — (tt, 0) 
varies for different cuprate families as expected from 
the above mentioned properties of the tt't" J model dis- 
persion and the values of t' and t" predicted by band 
theory This state of affairs, namely the agreement 
between tt't" J model predictions and experimental ob- 
servations, offers support to the relevance of this model 
in the underdoped regime of cuprates. 

2. Non-zero hole density 

The aforementioned experimental results concern 
ARPES measurements on undoped cuprate samples. 
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There is, however, experimental evidence that the dis- 
persion present in these samples persists as a broad high 
energy hump even away from half-filling and in the pres- 
ence of SC long-range order. 47 i 48 i 50 In Refs. lilHi this 
high energy dispersive feature is paralleled to the e^ k 
band in the "doped carrier" MF theory, which up to cor- 
rections of order x is given by the bare dopon dispersion 
7fc and, thus, by ti and t%. 

The above hump energy at (0, it) and (ir, 0), also known 
as the high energy pseudogap, lowers continuously as the 
hole concentration is increasedi 4 ^^ In order to reproduce 
this experimental evidence within the "doped carrier" 
MF approach the effective hopping parameters t% and t$ 
must be doping dependent. In particular, note that if 
<2 = 2<3 = tj the resulting dispersion is flat along the 
(0,7r) — (tt, 0) line. Therefore, in order to reproduce the 
decrease of the high energy pseudogap scale, the t% and t^ 
doping dependence must come from the doping induced 
renormalization of i' and t" . As a result, in the presence 
of non-zero hole density we use 

t 2 = tj + r(x)t' 

t 3 = ^+r(x)t" (35) 

instead of Expression i|34|l . Here, r(x) is a renormal- 
ization parameter which satisfies r(0) = 1 and that de- 
creases with increasing x. 

The above doping induced renormalization of t' and t" 
is suggested from comparison to experiments. However, 
below we argue that such a renormalization is consistent 
with theoretical studies of the tt't" J model. In Sec. IIII Bl 
we show that the dopon dispersion and, thus, e^~ fc as well, 
depends on the local spin correlations that enter the av- 
erage (|25fl . In the present MF approach this average is 
set by hand and is not calculated self-consistently. There- 
fore, it misses the doping induced changes in the underly- 
ing local spin correlations. In Ref. 0it is explicitly shown 
that, in the tt't" J model, the spin correlations induced 
by a hole hopping in a lattice of antiferromagnetically 
correlated spins strongly frustrate t' and t" . Hence, we 
expect that upon calculating the spin average l|25[l self- 
consistently the above renormalization of intra-sublattice 
hopping processes is properly reproduced. 

The just mentioned theoretical results indicate that 
the renormalization coefficient r(x) should decrease with 
x, however, they give no information toward its explicit 
functional dependence. We thus recur to experimental 
data which indicates that the high energy pseudogap 
scale vanishes around x « 0.30^^ to chose r(x) to van- 
ish at x = 0.30. In addition, we consider r(x) to interpo- 
late linearly between its x = and x = 0.30 values, which 
yields r(x) = (l — We have also considered alterna- 
tive interpolation schemes (not shown), say by changing 
the exponent of (l — ^) from 1 to 2, without affecting 
our general conclusions. 

We remark that, even though t 2 and t 3 are used to 
control the high energy dispersion t2,k hi consonance with 
experiments, there is no such direct experimental input 



on the low energy band e± t k, and all its prope rties result 
from the theory. Interestingly, Refs. I43u52ll53l find that a 
variety of low energy spectral properties associated with 
the ei,fe bands are consistent with ARPES and tunneling 
experiments on both hole and electron doped cuprates. 
We note that the cuprate hole doped (HD) regime can 
be addressed using t' « — 2t" « — J^i^ which within 
the context of the "doped carrier" MF theory reduces to 
using the effective hopping parameters 

f-'-^-is) 

In the electron doped (ED) regime t' and t" change 
sign 7 *— and, hence, in this case t^ and £3 become 

<f-.r + j(i-£) 
"f-i-^-is) ,37 » 

F. Staggered magnetization decoupling channel 

Both theoretical^ and experimentalSiSiZSiZi evidence 
support that the above MF theory Hamiltonian J3UJ), 
which assumes a spin liquid background, breaks down 
at and close to half- filling, where long-range AF order 
sets in. Therefore, here we extend the "doped carrier" 
MF approach in order to account for the staggered mag- 
netization decoupling channel. Specifically, we introduce 

m= (-iy* +i v(S?) (38) 

and 




u=2,3 u£uNN 

which are the lattice spin staggered magnetization and 
the dopon staggered magnetization respectively. 

The contribution from the above decoupling channels 
adds to l|3UI) so that we obtain the new MF Hamiltonian 
which allows for the presence of the AF phase 

H%f = H™/ + 2J*Nm 2 - ANmn- 

- 2 (J*m -n)J2 ^L^.tt)^- 

k 

- 2m^(7 fe + ^d)vl + ^)Vk (40) 

k 

It is well known that the above MF AF decoupling 
scheme overestimates the strength of magnetic moments. 
To effectively include the effect of fluctuations, which de- 
crease the staggered magnetization, we introduce a renor- 
malized exchange constant J* — XJ in the staggered 
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FIG. 1: The "doped carrier" MF phase diagram includes 
the states: AF, d-wave SC (dSC), strange metal (SM), Fermi 
liquid (FL), and pseudogap with and without Nernst signal, 
labeled by N and PG respectively. The hopping parameters 
(a) ti = t = 3 J and (b) ti = t = 1.5J are used, together with 
1361 and 1371 for the hole doped (HD) and electron doped 
(ED) regimes respectively. 



magnetization decoupling channel^ 2 - The renormaliza- 
tion factor is determined upon fitting the MF staggered 
magnetization at half-filling to the quantum Monte-Carlo 
estimate m = 0.31 iS In the present case, this condition 
requires A = 0.34i£i 

IV. DOPED CARRIER MEAN-FIELD PHASE 
DIAGRAM 

Starting from the MF Hamiltonian i jlDj l the MF phase 
diagram can be computed for different values of doping 
x and temperature T by requiring the self-consistency 
of the MF parameters and by determining the value of 
the Lagrange multipliers fid and ao that enforce the dop- 
ing density ((d\di) = x) and the global SU(2) projection 
((/»•/») — 1) constraints. If we ignore states with coexist- 
ing AF and SC orders we can separately consider those 
cases when both bo,bi ^ and those cases when both 
to, n 7^ 0. The corresponding saddle-point conditions 
can then be cast analytically, as shown in Appendix [B] 
for m,n — and in Appendix ICl for bo, b\ = 0. 

The generic "doped carrier" MF theory phase diagram 
was first computed in Ref. 13- Here, we show the MF 
phase diagram for a new set of parameter values, namely 
for NN hopping t x =t = 2>J [Fig. QJa)] and t\=t = 1.5 J 
[Fig. HJb)] and for the intra-sublattice hopping parame- 
ters ti and t3 in Expressions H36|l and H37JI. We remark 
that t = 3J together with (|36(l describe the parameter 
regime of relevance to HD cuprates and that t = 3 J to- 



gether with (|37|l address the ED regime. We consider the 
t = 1.5 J case in order to illustrate the role of t on the 
local energetics (see Sec. II V Ajl . 

Fig. shows that the MF phase diagram includes six 
different regions 7 ^ which correspond to distinct physical 
regimes. These regimes have been discussed within the 
context of slave-boson MF theorjiiSiSiZ 6 . and, in what 
follows, we briefly review their properties: 

(i) Antiferromagnet (AF) - (bo), (b\) — and m,n 7^ 
0. At and close to half-filling the lattice spins form local 
staggered moments and, thus, the low energy spin exci- 
tations are spin waves. Charge carriers move within each 
sublattice and form small Fermi pockets whose volume 
equals the doping level. 

(ii) Strange metal (SM) - (bo), A, m, n = 0. In 
this spin liquid state the low lying excitations are spinons, 
which have an ungaped Fermi surface and, as such, lead 
to a large low energy spin density of states. Since spinons, 
which are charge neutral, do not coherently mix with do- 
pons, which are charged, the resulting phase is an inco- 
herent metal. 

(Hi) Pseudogap metal (PG) - (\bo\), (\bi\), to, n = 
and A ^ 0. When spinons pair up in the d-wave channel 
a gap opens in the uniform susceptibility in agreement 
with the observed reduction of the Knight shift in the 
pseudogap regime. 77 Despite the gap, spin correlations 
at (tt, 7t) are enhanced by the gapless U(\) gauge field, 62 
as expected in the underdoped regime close to half-filling. 

(iv) Nernst regime (N) - (bo),(bi),m,n = and 
(\bo\), (\bi\), A 7^ 0. Below the MF spinon-dopon pairing 
temperature (|&o|)) ^ and the motion of spinons 
leads to a backflow in the charged &o,i fields so that, effec- 
tively, spinons transport electric charge. Consequently, 
in the presence of d-wave spinon pairing (A ^ 0) the 
system displays cZ-wave SC correlations. Since the mag- 
nitude of the spinon-dopon pairing field vanishes toward 
half-filling, in the underdoped regime phase fluctuations 
may prevent the onset of true long-range SC order, lead- 
ing to a region in phase space where (bo), (b\) = even 
though (|6o|), (\bi |) 7^ 0. In this region, which we call the 
Nernst region, short-range SC fluctuations can be exper- 
imentally detected through the Nernst effect 7 ^ 7 ^*^^ or 
the diamagnetic response^ 

(v) d-wave superconductor (dSC) - m,n — and 
(&o),(&i),A 7^ 0. Below the Kosterlitz-Thouless tran- 
sition temperature for the above mentioned phase fluc- 
tuations, the fields bo and b\ display long-range phase 
coherence and the electronic system is a d-wave super- 
conductor. 

(vi) Fermi liquid (FL) - A, to, n = and (bo), (bi) =/= 0. 
Since spinons, which are not superfluid (A = 0), hy- 
bridize with dopons they effectively are charged spin-1/2 
fermionic excitations with a large Fermi surface. There- 
fore, the electronic system is in the Fermi liquid state, as 
expected from the Ioffe-Larkin sum rule^ 

Fig. ^a) depicts the MF phase diagram in the param- 
eter regime of interest to the cuprates and shows that 
antiferromagnetism is very feeble on the HD side leaving 
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FIG. 2: Plot of the doping level x a f at which d-wave SC 
order replaces AF order (circles) and of the maximum doping 
level Xsc of the d-wave SC dome (crosses) as a function of 
t' = -2t" = -aJ for both t = 3.0J (solid line) and t = 1.5J 
(dash-dot line). Note that a — 1 and a = — 1 reproduce 
the hole doped parameters in l|36jl and the electron doped 
parameters in 137H respectively. 



room for a large pseudogap region, which is mostly cov- 
ered by the Nernst regime^i and for a large SC dome. 
In the ED case, however, AF order is more robust and 
covers a considerable fraction of the SC dome, which is 
smaller than for the HD regime, as well as nearly all 
of the Nernst region, in conformity with the lack of ex- 
perimental evidence for a vortex induced Nernst signal 
in these materials^ In fact, the above MF phase dia- 
gram not only captures the asymmetry between the HD 
and ED regimes, in agreement with previous numeri- 
cal studies^ii but is semi-quantitatively consistent with 
that of real materials. Indeed, based on the experimental 
and numerical input referred to in Sees. IIII El and IIII Fl 
we find that: the SC dome extends up to x w 0.24 on the 
HD side while coming to an end at x « 0.16 for the ED 
regime; the maximum T c is ~ 100K in the HD case while 
it is smaller on the ED side. Furthermore, the renormal- 
ized J* , which controls the strength of local magnetic mo- 
ments, is determined by the behavior of the MF theory 
at half-filling and, yet, it correctly predicts that AF order 
ceases to exist at a doping level which is consistent with 
experiments on both HD and ED compounds £§2/22/2^ 

A. The role of t, t' and t" 

The phase diagrams in Fig. ^ show that as we move 
away from half-filling the AF phase is initially replaced 
by a state where spins are paired in the d-wave singlet 
channel. This d-wave gapped spin liquid state enhances 
the doped carrier kinetic energy while preserving much 
of the magnetic exchange energy^ As the doping level 
increases carrier motion further frustrates the exchange 
energy and eventually closes the d-wave gap. Therefore, 
the doping evolution of phases is closely connected to the 
doped carrier dynamics and, below, we explore the micro- 
scopic picture that underlies how the hopping parameters 



t, t' and t" affect the robustness of AF and SC correla- 
tions (see the different phase diagrams in Fig. ^as well 
as how x a f and x sc , which stand for the doping level at 
which d-wave SC order replaces AF order and the maxi- 
mum doping level of the d-wave SC dome, change with t 
and t' = -2t" in Fig. |2J). 

The combined effect of the t, t' and t" terms may en- 
hance or deplete the hopping between first, second and 
third NN sites and, thus, increase or decrease the cou- 
pling between doped carriers and specific surrounding 
spin correlations. The particular case of NN hopping is 
controlled by t and it frustrates the spin correlations in- 
duced by the exchange term in the Hamiltonian. Hence, 
increasing t/J enfeebles antiferromagnetism and reduces 
the d-wave spin pairing amplitude A, as supported by 
the decrease of both x a f and x sc in Fig. [2] when we go 
from t = 3 J to t — 1.5J. 

The hopping between second NN sites is controlled 
both by t and t - if the latter parameter is positive then 
t and t' processes interfere constructively to enhance sec- 
ond NN hopping, while if t' < these processes inter- 
fere destructively to deplete second NN hopping. The 
same argument applies to third NN hopping if we take 
t" instead of tf above. Since both second and third NN 
hopping do not harm the staggered spin configuration of 
the AF state, increasing either t' (as in the ED side) or 
t" (as in the HD side) stabilizes AF order and increases 
x a f (Fig. |2J) ■ t' in the ED regime is larger than t" in 
the HD regime and, therefore, this effect is more promi- 
nent in the former case, leading to the aforementioned 
asymmetry in the phase diagrams. 

Figs, n an d 121 support that the values of t' and t" 
also affect how doped carrier motion couples to d-wave 
SC correlations. To understand this effect note that, in 
a d-wave superconductor, the condensate induces intra- 
sublattice hopping processes where the amplitude for a 
hole to hop between second NN sites is negative while 
the amplitude for a hole to hop between third NN sites is 
positive. These processes, and thus SC order as well, are 
frustrated when t' — —2t" > 0. A different way to picture 
the above argument is to note that when t' = —2t" > 
the doped carrier dispersion is gaped at (§ , f ) , in which 
case, a vacancy hopping in the presence of local AF 
correlations frustrates the spinon d x 2_ y 2-w&ve dispersion 
and weakens superconductivity. Therefore, x sc decreases 
when t' and t" vary between the HD and ED regimes 
(Fig. [21 • The same applies to the highest SC T c , in 
agreement with experiments and other theoretical ap- 
proaches. Indeed, band theory calculations together with 
experimental data support that the maximum T c for var- 
ious cuprate families increases with —t'/tA ARPES re- 
sults also suggest the correlation between the high en- 
ergy pseudogap scale, which is controlled by t' and t" , 
and the maximum T C J& In addition, variational Monte 
Carlo calculations further substantiate the above role of 
t' in determining the robustness of SC correlations^ 

The previous digression on the roles played by t, t' 
and t" in the interplay between spin and charge degrees 



12 



of freedom may be of interest to understand the striking 
differences between as-grown and oxygen reduced elec- 
tron doped cuprate compounds. The former samples are 
not SC and display long-range AF order up to x ~ 0.20. 86 
After the oxygen reduction process, which removes about 
1% of the oxygen atoms in these materials, AF order is 
destroyed at x ~ 0.10 — 0.14 and, at higher doping val- 
ues, the samples superconducti22i£i*Sl Within the present 
context, we propose this sharp change follows the alter- 
ation of the effective in-plane parameters t, t' and if'. 
Since both if and t" depend on the chemical composi- 
tion outside the copper-oxide layers^ if such a composi- 
tion changes in a way that the magnitude of if and t" 
decreases, the phase-space volume of the AF phase is re- 
duced while SC correlations are enhanced. Alternatively, 
if the oxygen reduction process acts on the copper-oxide 
planes in such a manner that t effectively decreases, then 
4 ~ y increases, which favors superconductivity over an- 
tiferromagnetism. The latter scenario receives support 
from experimental evidence for the removal of oxygen 
atoms from the copper-oxide planes under the reduction 
process in both PCCOS and NCCOS 



V. SUMMARY 

In this paper we explicitly derive a recently 
introduced^ formulation of the tt't" J model in terms of 
projected dopon and spin operators instead of projected 
electron and spin operators. Since dopons describe the 
carriers doped in the half-filled system, we name it the 
"doped carrier" formulation of the tt't" J model. Close 
to half-filling the doped carrier density is small and, thus, 
we circumvent the "no-double-occupancy" constraint. In 
particular, we propose that the effect of the projection 
operators V in the usual formulation of the tt't" J model 
l[TJl is captured by the interaction between doped carri- 
ers and lattice spins in \Yi\. This interaction explicitly 
accounts for the interplay between local spin correlations 
and the hole dynamics in doped Mott insulators. 

The tt't" J model Hamiltonian in the enlarged Hilbert 
space [Expression l|10f> ] provides a new starting point to 
deal with doped spin models, which we pursue to de- 
velop a new, fully fermionic, MF theory of doped Mott 
insulators. The resulting "doped carrier" MF theory is 
constructed to address the low doping and low temper- 
ature regime of the tt't" J model, and properly accounts 
for the frustration of NN hopping due to the strong lo- 
cal AF correlations present in such a limit. Since a hole 
hopping in an antiferromagneticaly correlated spin back- 
ground induces new local spin correlations that strongly 
renormalize t' and t" ^ in (|35|l we introduce a phe- 
nomenological doping dependent renormalization factor 
r(x). Specifically, we choose r(x) so that the high en- 
ergy MF dispersion k reproduces the evolution of the 
high energy pseudogap scale at (0, ir) observed in ARPES 
experiments!^^ Remarkably, using t, if and t" moti- 
vated by band theory^S& and with the aforementioned 



little experimental input the "doped carrier" MF theory 
leads to a semi- quantitative correct phase diagram for 
both HD and ED cuprates [Fig. Ha)]. In particular, 
in the HD case a large SC dome and extended pseudo- 
gap regime are obtained, while superconductivity is much 
weaker on the ED side where it is partly overtaken by the 
robust AF phase. 

In the low doping limit of the generalized-i J model 
it is meaningful, and we believe useful as well, to think 
of vacancies encircled by local moments. The "doped 
carrier" approach then provides a framework to ad- 
dress the hole dynamics in the presence of various lo- 
cal spin correlations. In the hereby developed MF the- 
ory, an effective multi-band description captures the 
effect of two different local correlations, namely stag- 
gered moment and d-wave singlet bond correlations: the 
vacancy in the one-dopon state is encircled by stag- 
gered local moments which inhibit NN hopping; upon 
spinon-dopon mixing the vacancy changes the surround- 
ing spin background and gains kinetic energy, which is 
the driving force for the d-wave SC state at low dop- 
ing. The above two-band description receives support 
from more rigorous calculations on the two-dimensional 
tJ and Hubbard models) 8 i 9 i 20 i 21 i 22 i 23 and is consistent 
with the various spectral features observed by ARPES 
experimentsi^*2iL&2i£i We remark that in the present ap- 
proach superconductivity arises due to the change im- 
posed on spin correlations by the motion of doped carri- 
ers. This process is captured by the spinon-dopon mix- 
ing, which can be viewed as a new mechanism for super- 
conductivity. We believe this mechanism to be relevant 
for the experimentally obtained low values of J/t. In the 
large J/t regime a distinct mechanism, namely the min- 
imization of the number of missing AF links, may lead 
doped carriers to form local pairs, as proposed by the 
"antiferromagnetic-van Hove model" 

Since the "doped carrier" approach addresses the in- 
terplay between the doped carrier dynamics and different 
background correlations, we can discuss how the strength 
of the latter depends on the hopping parameters. Inter- 
estingly, we find that short-range AF correlations can 
enhance or deplete d-wave SC correlations depending on 
the sign oit' and t" . Specifically, if these hopping param- 
eters favor a gap in the single hole tt't" J model dispersion 
at (^, ^), as in the ED regime, the vacancies when in the 
presence of local AF configurations frustrate the <i-wave 
SC gap and, thus, superconductivity as well. The oppo- 
site effect occurs if, instead, if and t" induce a gap at 
(0, 7r) in the single hole tt't" J model dispersion. Hence, 
the interplay between short-range AF and <i-wave SC cor- 
relations provides a microscopic rational for the role of 
t' and t" in strengthening superconductivity as expected 
by experiments and other theoretical approaches ££L§k 
In this context, note that the "doped carrier" MF theory 
not only accounts for the hole/electron doped asymme- 
try but also offers possible scenarios for the difference 
between phase diagrams of as-grown and oxygen reduced 
electron doped samples 
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The "doped carrier" formulation provides a MF the- 
ory to describe superconductors with strong local AF 
correlations due to a proximate Mott insulating state. 
The signature of such correlations is explicit in a vari- 
ety of cuprate experimental data that deviates from the 
pure BCS behavior. In particular, Refs. I43TI52TI531 show 
that several non-trivial features of the electron spectral 
function and of the tunneling conductance spectrum of 
cuprates are reproduced by the herein presented MF ap- 
proach. 
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APPENDIX A: RELATION TO SLAVE-BOSON 
APPROACH 

In the present paper a new formulation of the tt't" J 
model is used as a starting point to develop a new MF 
theory of doped Mott insulators. This "doped carrier" 
MF theory, which to the authors' best knowledge dif- 
fers from other MF theories in the literature, bears some 
relations to the slave-boson MF approach to the same 
modeL 90 Both approximations recur to the fermionic 
representation of spin operators and, thus, in the un- 
doped limit both MF approaches are equivalent. How- 
ever, these MF theories deal with doped carriers in dis- 
tinct manners since they introduce different operators 
to account for charged degrees of freedom in doped sys- 
tems - the slave-boson formulation introduces the spin- 
less charged bosonic holon operator h\ = [/1J1/1J2] while 
the "doped carrier" formulation introduces the spin- 1/2 
charged fermionic dopon operator d\. In this appendix 
we clarify the relation between holon and dopon opera- 
tors and, thus, the relation between the extensively used 
slave-boson approach and the new "doped carrier" for- 
malism. 

Within slave-boson theory, the projected electron op- 
erators, which are the building blocks of all other physical 
operators, can be written in terms of holons and spinons 
as 



t <D_ 



T2 



(Al) 



as long as we constrain ourselves to the physical Hilbert 
space defined by lM(ripi+h\trhi) = O^MThis formula- 
tion introduces an SU(2) gauge structure since physical 



operators are invariant under the local transformation 
tpi — > Witpi and h{ — > Wihi, where Wi is any SU(2) ma- 
trix. Hence, both spinons and holons carry an SU (2) 
gauge charge in addition to the physical spin and electric 
charge quantum numbers. 

Interestingly, in the "doped carrier" formalism dopons 
and spinons can form singlet pairs, as captured by the 
operators B^j introduced in Sec. IIII Bl These are elec- 
trically charged spinless bosonic fields which also carry 
the above SU (2) gauge charge and, therefore, they have 
the same quantum numbers as holon operators in the 
slave-boson framework. To make the relation between 
holons and spinon-dopon pairs explicit it is convenient 
to rewrite the projected electron operators in terms of 

dopons and spinons. We use Sf = | (f}^fi,l — fjifi,l), 

5 2 + = /j t/i.j. and — f}ifi,i to express the operators 
in as 



lY/r = —//r (f iA d itl - f iA d iA )v 



Vc\{P = -j=]\ A V (f iA di,i - f iA d iA )V (A2) 



Since the above equality only holds for fjfi 
be recast as 



1 it can 



Vc h V 
Vc li V 



1 



f .r{.r; ,i,.- ■ //. </,. v 1 



\hAV (/J^.t + /;.//,..) V (A3) 



V2 



Using the slave-boson formulation language, in the phys- 
ical Hilbert space there can be at most one holon per 
site. Say it happens to be a hi. 1 holon. Then, the pro- 
jected electron operators in (|A1|) resemble those in (|A2|) 
with the hi t i holon replaced by the spinon-dopon pair 
— fi,[di,i). If, instead, there exists a holon 
on site i, the projected electron operators in IjAlfl re- 
semble those in (|A3|I with the hi t 2 holon replaced by the 

spinon-dopon pair yfj^di^ + fj^d^ij . In this sense, the 

hi t i and hi 2 holon operators are related to specific singlet 
pairs of spinons and dopons. 

The above correspondence between holons and spinon- 
dopon pairs can be used to compare MF phases in the 
slave-boson and "doped carrier" approaches. Specifically, 
for a given spinon state, the physical symmetries of a 
phase described by a certain pattern of holon condensa- 
tion (h\) = [{h\i){h\ 2 )] in the slave-boson formulation 
are the same as those of a state where spinon-dopon pair- 
ing in the "doped carrier" formulation yields 



Bio 



-(h i<2 ) 
(4i> 



{hi,: 
{hi 



(A4) 



This result shows that holon condensation in the slave- 
boson formalism leads to the same phases as the spinon- 
dopon pairing transition observed in the "doped carrier" 
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approach. In particular, in the presence of d-wave paired 
spinons described by ansatz Q2U|). the rf-wave SC state 
obtained upon the condensation of holons = [0 ho] 
is the same as the one that occurs in the presence of the 
spinon-dopon pairing described in l|29|) . 

To understand the connection between holons and 
spinon-dopon pairs it is useful to consider the local pic- 
ture of these objects. The dopon is an entity which car- 
ries charge and spin and is formed by the vacancy plus 
a neighboring spin. More precisely the dopon spin is 
carried by the staggered local moments that surround 
the vacancy in the one-dopon state. Since these local 
moments frustrate doped carrier hopping the spin state 
encircling the vacancy in the one-dopon state is mod- 
ified to optimize the doped carrier kinetic energy. In 
the "doped carrier" MF theory this interaction between 
doped carriers and the surrounding spins is captured by 
the ~ <ft ffd term in (1221) which can be recast as ~ bf'd. 
This term drives the decay of the dopon into a spinless 
spinon-dopon pair and a chargeless spinon. Physically, 
this process means that the spin background forms a sin- 
glet with the dopon spin thus altering the spin configura- 
tion around the vacancy. The vacancy is then encircled 
by a local spin singlet configuration, which corresponds 
to the local picture of a holon. The doped carrier spin- 
1/2 is absorbed by the spin background in the form of a 
spinon excitation. 

The above picture suggests that the holon is a compos- 
ite object and that the "doped carrier" approach captures 
its internal structure. Such an internal structure should 



then be apparent in the electronic spectral properties. In 
Ref. these properties are discussed at length and the 
holon internal structure is argued to be reflected in the 
momentum space anisotropy (also known as the nodal- 
antinodal dichotomy^) of the electron spectral function. 



APPENDIX B: PARAMAGNETIC MEAN-FIELD 
SELF-CONSISTENCY EQUATIONS 

Below we write the set of self-consistent conditions for 
the paramagnetic MF Hamiltonian l|30|l which are deter- 
mined by the saddle point equations 



dF 
dx 



dF 
OA 



dF_ 
db~o 



dF_ 

db[ 



dF 



dF 







(Bl) 



where F is the MF free-energy of the fermionic system 
described by (|3Ufl . namely 



|^^( x 2 + A 2 )-2&o&i-Ml-*)- 



N ^ 

k 



1 + cosh 



T 



1 + cosh 



c 2,fe 



T 



(B2) 



The resulting self-consistent equations for the mean- 
fields x, A, bo and b±, as well as the doping concentration 
and SU(2) projection constraints, are: 



1 _ i 2^/32 + a£ 2/3 2 + Kr + Kr- 7 2 
x = -J^l^( cosk x + cosk y) \ a l A k 



2V^fc 



Bi. 



(B3) 



A = — 



J^J2 a k ( cosk x - cos k y ) \ A k 

k I 



2/3g + K) 2 + K) 2 - 7 I 



Bi 



(B4) 



&i = — — - /3fe [ti (cos k x + cos k y ) + 2t 2 cos k x cosk y + £3 (cos 2^ + cos2k y )] 



f Afe _(7 fe + a|f + K) 2 5fe ' 
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k K. 
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k 



B k 



(B8) 



where we introduce 



and the self-consistency equations 



sinh | 



and 



-i,fe 



1 + cosh | 



sinh | 



-i,fe 



1 + cosh | ^ 



sinh | -|£ 



c 2 : fc 



1 + cosh [ ^ 
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sinh | -|£ 



c 2,fe 



1 + cosh 



~2,fc 



The remaining notation is defined in l|31|) and l|33[l . 



(BIO) 



reduce to 



dF _dF _dF _dF _ dF 
dx dA dm dn d[id 



(C3) 



2N 4^ 



sinh [ s ' k 
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1 + cosh 
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APPENDIX C: ANTIFERROMAGNETIC 
MEAN-FIELD SELF-CONSITENCY EQUATIONS 



In Sec. IIII Fl we extend the "doped carrier" MF theory 
to include the staggered magnetization decoupling chan- 
nel and obtain the MF Hamiltonian l|40|) . In this paper 
we do not consider states with coexisting AF and SC or- 
der and, in what follows, we set bo and b\ to zero in (|4Up . 
Since in this case spinons and dopons do not mix we can 
define two spinon and two dopon bands, namely 



= fc = ±v /K) 2 + K) 2 + ^ 

e d,k = (1 =F 2jm|) 71, — A*d 



(CI) 
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where isp. 



-2(J*m-n). 



In the absence of 



spinon-dopon hybridization ag 
- (fx - ^r) (cosk x + cosk y ). 
The resulting MF free-energy is 



and a 
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